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ABSTRACT 

We present Hubble Space Telescope observations of the upper part {Tes > lO"* K) of the white 
dwarf cooling sequence in the globular cluster 47 Tucanae and measure a luminosity function of hot 
white dwarfs. Comparison with previous determinations from large scale field surveys indicates that 
the previously determined plateau at high effective temperatures is likely a selection effect, as no such 
feature is seen in this sample. Comparison with theoretical models suggests that the current estimates 
of white dwarf neutrino emission (primarily by the plasmon channel) are accurate, and variations are 
restricted to no more than a factor of two globally, at 95% confidence. We use these constraints to 
place limits on various proposed exotic emission mechanisms, including a non-zero neutrino magnetic 
moment, formation of axions, and emission of Kaluza-Klein modes into extra dimensions. 

Subject headings: astroparticle physics - dense matter - elementary particles - neutrinos - stars: 
luminosity function - white dwarfs 


1. INTRODUCTION 

Despite its modest beginnings in a sheepish proposal 
by Pauli in his famous letter of 19310, the neutrino has 
emerged as an object of great interest to the contempo¬ 
rary particle physics community. The interpretation of 
solar and atmospheric neutrino data (Davis, Harmer & 
Hoffman 1968; Fukuda et al. 1998) in terms of neutrino 
oscillations requires that neutrinos have finite mass and 
exhibit flavour mixing, which is one of the few pieces of 
evidence that currently points to the existence of phys¬ 
ical phenomena beyond the standard model of particle 
physics. However, the small interaction cross-section of 
neutrinos means that they are quite difficult to produce 
and to study in terrestrial experiments. On the other 
hand, certain classes of stars produce neutrinos in copi¬ 
ous quantities and neutrino cooling plays an important 
role in certain aspects of stellar evolution. 

Indeed, when central densities p > 10'^g.cm~^ and cen¬ 
tral temperatures are ~ 10® K, the stellar luminosity in 
neutrinos can exceed that in photons (Fowler & Hoyle 
1964; Vila 1966; Lamb & Van Horn 1975). Such condi¬ 
tions are realised during the core evolution of intermedi¬ 
ate mass stars, and neutrino cooling regulates the rate of 
evolution of stars as they transition from the Asymptotic 
Giant Branch to the top of the white dwarf cooling track. 
Models for masses ~ 0.53Mq, an appropriate mass for 
those white dwarfs forming in old globular clusters today, 
show that for the first ^ 5 x 10® years of the white dwarf 
existence, the rate of cooling is dominated by neutrino 
emission. Thus, if we can amass a large enough sample 
of white dwarfs with ages <10 Myr, the study of the re¬ 
sulting hot white dwarf luminosity function can provide 
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a direct test of our understanding of neutrino emission 
in astrophysical settings. 

Such studies became a possibility with the construction 
of a luminosity function for hot white dwarfs in the solar 
neighbourhood from the Palomar-Green survey (Flem¬ 
ing, Liebert & Green 1986). The resulting comparison 
with models for the population evolution of young white 
dwarfs enabled tests of our understanding of neutrino 
emission, including constraints on the possibility of extra 
cooling from hot cores due to the existence of processes 
beyond the standard model, such as axion emission or the 
existence of a neutrino magnetic moment (Wang 1992; 
Blinnikov & Dunina-Barkovskaya 1994). However, the 
short evolutionary timescales and relative rarity of hot 
white dwarfs meant that these constraints were not im¬ 
proved until the advent of the modern era of large scale 
sky surveys, such as the Sloan Digital Sky Survey (SDSS 

- York et al. 2000) or SuperGosmos Sky Survey (SSS 

- Hambly et al. 2001). This has led to a resurgence of 
interest in this question, with improved constraints on 
the hot white dwarf luminosity function presented by De 
Gennaro et al. (2008), Krzesinski et al. (2009) and Row¬ 
ell & Hambly (2011) and resulting theoretical analyses 
(Isern et al. 2008; Miller Bertolami 2014). 

However, the construction of a luminosity function 
from an all sky survey must confront a variety of selection 
effects, including a non-uniform selection for spectro¬ 
scopic follow-up and biases in the assignment of absolute 
magnitudes due to distance uncertainties (see Krzesinski 
et al. 2009 for a discussion). As a result, we describe here 
the construction of a hot white dwarf luminosity function 
from an alternative source, namely the globular cluster 
47 Tucanae. This project is a follow-up to a previous 
deep observation to constrain the luminosity function of 
cool white dwarfs (Hansen et al. 2013; Goldsbury et al. 
2012) but with a wider field coverage (albeit to shallower 
depth) and UV exposures, in order to increase the num¬ 
ber of bright white dwarfs detected. This program has 
an advantage over all-sky surveys in that all the white 
dwarfs lie at the same distance and are selected in a uni- 


2 


Hansen et al 


form manner, albeit with its own selection effects, to be 
discussed below. As such, our resulting luminosity func¬ 
tion provides a constraint on hot white dwarf physics 
that is completely independent of that derived from the 
stars in the solar neighbourhood. 

We describe below (§12]) two sets of observations taken 
with the Hubble Space Telescope in 47 Tucanae. In § [3] 
we compare these to a set of white dwarf atmospheric 
models and describe resulting construction of the hot 
white dwarf bolometric luminosity function. We also 
describe the construction of a set of theoretical cooling 
models (§11 used to explore the possible ranges of phys¬ 
ical inputs for comparison to the observations, including 
both explorations of the neutrino physics, possible non¬ 
standard contributions to the luminosity, and potential 
sources of confusion, such as uncertainties in the thick¬ 
ness of the Hydrogen layer on the surface of white dwarfs. 
Finally, in § [S] we compare the models and the data to 
place constraints on the various parameters discussed. 

2. OBSERVATIONAL DESIGN AND DATA 
ANALYSIS 

The observational design of this program (GO-12971; 
PI H. Richer) involves parallel imaging observations with 
HST’s WFC3 and ACS cameras. The data were collected 
from Nov 14th 2012 to Sep 20th 2013. The 10 orbits were 
split into 10 visits, each one with an orientation offset to 
ensure a 360 degree mapping of the cluster central region 
when all data are combined together. In this arrange¬ 
ment, the center of 47 Tuc was placed near the corner of 
all WFC3 observations, so the contiguous region from the 
mosaic extends approximately one WFC3 diagonal field 
in radial extent. The parallel ACS pointings approxi¬ 
mately abut the WFC3 helds and therefore complete an 
annulus around the cluster center. Taken together, this 
arrangement provides high-resolution HST imaging cov¬ 
ering the inner 8.5 arcmin of 47 Tuc. 

The data analysis was performed using the prescrip¬ 
tion and procedures described in Kalirai et al. (2012). 
Within each visit, the observations contain two WFC3 
exposures in F225W (380s and 700s) and F3361F 
(485s and 720s), and two ACS exposures in F435kF 
(290s and 690s) and F555VF (360s and 660s). For 
the WFC3 observations, we retrieved the _raw and _flt 
frames and corrected the data for CTE losses using the 
WFC3 pixel-based empirical CTE correction software 
(http: //www.stsci.edu/hst/wfc3/tools/cte.tools). Eor 
the ACS data, we retrieved the _flc (i.e., CTE corrected) 
frames directly from the MAST archive. For each fil¬ 
ter on the WFC3 data, we executed a first pass of Mul- 
tiDrizzle (Fruchter & Hook 2002) to generate geometric 
distortion corrected images. We then performed astrom¬ 
etry on the stars in these images (using the DAOPHOT H 
software; Stetson 1987; 1994) and calculated transforma¬ 
tions that mapped all images on the 10 visits to the same 
frame of reference. These offsets and rotations were then 
supplied to MultiDrizzle in a second pass as a “shift” 
file and a stacked (i.e., drizzled) image was produced in 
each filter for the entire 10 visit mosaic. Given the small 
overlap of the 10 pairs of images, we did not rescale the 
images or adjust the “pixfrac” parameter in MultiDrizzle. 
As the ACS data contains negligible overlap between the 
ten fields, each of the fields were treated independently 
to create ten drizzled images from the pairs. 
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Fig. 1.— The panel shows the WFC3 UV data on the core field, 
with the most prominent being the cluster main sequence and giant 
branch. A prominent white dwarf cooling sequence is also evident, 
and separated from the background SMC population by the dashed 
curve. The faint population of detections blueward of the white 
dwarf sequence are spurious. 

The single drizzled WFC3 image in each filter, and 
the ten drizzled ACS images in each filter, were sub¬ 
jected to PSF-fitted photometry and morphology using 
DAOPHOT H. As described in Kalirai et al. (2012), our 
procedure involved finding all sources that are at least 
2.5 ct above the local sky, selecting candidate PSF stars 
based on their brightness and isolation, generating a spa¬ 
tially variable PSF from these stars, and applying the 
PSF to all sources in the image using ALLSTAR. The 
resulting photometry from each filter was zero pointed 
to the VegaMAG system, and the photometry from the 
two filters on each camera was matched into a single cat¬ 
alog. 

2.1. The Core Field 

The CMD of all stars in the WFC3/UVIS drizzled im¬ 
age is shown in Figure [U and shows a clear white dwarf 
cooling sequence extending to F2251F ^ 25. At these 
short wavelengths, the main sequence of the background 
SMC is less prominent than at other wavelengths (e.g. 
Hansen et al. 2013; Kalirai et al. 2013; and ? 12.21) . The 
dashed curve shown in Figure [T] is used to identify the 
white dwarf population for the purposes of constructing 
a luminosity function. 

The approximate transition temperature between 
neutrino-dominated cooling and photon-dominated cool¬ 
ing is ~ 25,000K. A white dwarf with this temperature 
and log 5 = 8, at a distance of hq = 13.3 and a redden¬ 
ing of E(B-V)=0.04 (nominal 47 Tuc parameters), will 
have ESSOkF ~ 22.1. Our luminosity function contains 
216 white dwarfs brighter than this threshold. Although 
the exact value of this transition will vary depending on 
the distance, extinction and white dwarf parameters, it 
is clear that our data contains a substantial population 
of white dwarfs whose principal cooling mode is via the 
emission of neutrinos. 

To assess the photometric and astrometric error distri¬ 
bution of the final photometry in our images, as well as 
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Fig. 2.— This figure shows the data from the ACS visible CMD 
compiled from 9 separate pointings in the outer parts of the clus¬ 
ter. The three prominent populations are, from the right/red, the 
cluster main sequence, the SMC main sequence and the cluster 
white dwarf population. The dashed line indicates the CMD cut 
used to isolate the white dwarf population. We see that this sep¬ 
aration becomes more difficult at fainter magnitudes because the 
photometric scatter of the white dwarfs and the SMC background 
begin to overlap. 

the completeness of the data reductions, we performed 
artificial star tests. The artificial stars are modeled from 
the stellar PSF and scaled to reproduce the complete 
luminosity range of real stars in the drizzled images of 
each filter. The input stars were placed along each of 
the white dwarf and main sequence of 47 Tuc, and the 
input positions were set such that no two artificial stars 
overlapped each other on the image. As the mosaicing 
of our data set results in deeper observations near the 
center of the cluster (i.e., exposures from all visits over¬ 
lap in the core), we placed the artificial star tests over 
the complete radial range of our fields. The procedure 
to do this involved placing stars with the same intensity 
in radial grids. A total of 25 million artificial stars were 
used. The images with artificial stars were subjected to 
the photometric routines that were applied to the actual 
drizzled images, using identical criteria. The stars were 
recovered blindly and automatically cross-matched to the 
input starlists containing actual positions and fluxes. 

2.2. The Outer Fields 

Figure [2] shows the CMD from the ACS parallel fields. 
The presence of the SMC main sequence is much more 
prominent in these redder filters, and the photometric 
scatter is larger because the exposure time per object is 

10% that of the stars in the core data. As a result, the 
definition of the white dwarf cooling sequence becomes 
more ambiguous for F'bh'bW > 24.5, as the white dwarfs 
and SMC main sequence start to overlap. Nevertheless, 
we can still see a well sampled luminosity function for 
bright white dwarfs. 

The increased field coverage of the ACS data partially 
compensates for the lower stellar density relative to the 
core. The stellar density at ^ 3-4 arcminutes (the dis¬ 
tance of the parallel fields from the cluster center and 


slightly larger than the cluster half-mass radius - Trager, 
Djorgovski & King 1993) is more than an order of mag¬ 
nitude lower than that averaged over the WFC3 field 
(which encompasses ^2.7 arcmin), so that the increase 
in areal coverage by a factor of 10 results in the detection 
of 89 white dwarfs with F555VF < 23.5 (the magnitude 
for white dwarfs with Teff = 25, OOOK and the same dis¬ 
tance and extinction as used for the WFC3 estimate)^ 

We follow the same artificial star procedures as in S I2.1I 
to estimate completeness and photometric scatter. The 
radial dependance is much less in the ACS fields, because 
they lie at several core radii from the cluster, where the 
density gradient is much smaller across the field. 

3. FITTING ATMOSPHERIC MODELS TO THE 
DATA 

Our two data sets represent two observations of the 
same class of stars in multiple bandpasses and spatial lo¬ 
cations in the cluster. The underlying atmospheric mod¬ 
els should be the same in both bandpasses, and so we can 
hope to constrain the distance and extinction to the clus¬ 
ter based on fitting the position of the cooling sequence 
in both colour-magnitude diagrams simultaneously. 

In § [5] we will fit the data to models that include cool¬ 
ing histories, but we can compare the model and observed 
colours independently of the cooling history. This is be¬ 
cause the position of the cooling sequence in the colour- 
magnitude diagram depends on the model atmosphere, 
while the distribution of stars along the cooling sequence 
is determined by the rate of cooling. To constrain the 
former alone, we calculate the mean white dwarf colour 
as a function of magnitude (in 0.25 magnitude bins) in 
both the core and swathe helds, and fit an atmosphere 
model for a range of assumed white dwarf mass, extinc¬ 
tion and distance, using the artificial star tests to es¬ 
timate photometric scatter as a function of magnitude. 
We also considered various values of the surface Hydro¬ 
gen layer mass, which can inflate the radius of hot white 
dwarfs. Eigures [3] and |4] shows the resulting best fits in 
both cases. 

For the ACS data, the best ht distance is iiq = 13.21 ± 
0.16, for an extinction E{B — V) = 0.04 ± 0.02. The 
best fits are obtained for masses ^ 0.51 — O.53M0, and 
Hydrogen envelope mass fractions qn = 10“^. The er¬ 
ror bar on this determination is larger than our previous 
measurements in Woodley et al (2012) and Hansen et 
al. (2013) because the more limited magnitude range of 
the data allows for more covariance between distance, 
extinction and mass than in those analyses. The value 
here is slightly lower than in those cases 13.36 ±0.08 and 
13.32 ±0.09 respectively), but still consistent within the 
error bars. The constraint will be improved in § [5] when 
we use the full information in the cooling models. 

The best ht distance for the core data is 13.08 ± 0.08, 
again favouring white dwarf masses ^ 0.51 — 0.53Mq but 
now the best hts are obtained for low extinctions and 
large Hydrogen layer masses. Indeed, inspection of the 
comparison between models and data in Figure [4] sug¬ 
gests that the atmosphere models are not a good repre¬ 
sentation of the data in these passbands, as they are too 
blue at magnitudes F336VF > 23. Furthermore, the hts 
favour low extinction because the colour asymptotes to 
a value at bright magnitudes that is a poor match to the 
observed values, which drags the htted extinction down 
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Fig. 3.— The small open points represent the observed white 
dwarf population in the swathe of outer fields observed with the 
ACS camera. The large red points show the best fit atmosphere 
model, for a white dwarf of mass 0.53Mq, distance modulus po = 
13.21 and E(B-V)=0.04. The best fit is determined by minimising 
the difference between the model colour and the mean observed 
white dwarf colour at that magnitude. 



Fig. 4. — The small open points represent the observed white 
dwarf population in the core of 47 Tuc, observed with the WFC3 
camera. The large red points show the atmosphere model, for a 
white dwarf of mass O.SIMq, distance modulus po = 13.2 and 
E(B-V)=0.02. These are chosen to represent the most plausible fit 
consistent with other determinations of distance. We see that the 
models fit the data well for F336W< 23 but are too blue relative 
to the data for fainter magnitudes. It is this mismatch that drives 
the best fits to lower distances. 


and increases the mismatch at the faint end. The mod¬ 
elling of the core white dwarf sequence is further com¬ 
plicated by the effects of spatially varying completeness 
and the possibility of radial diffusion of stars relative to 
dynamical relaxation. However, neither of these effects 
can explain the discrepancy here, because the best-fit 
colour is determined at each magnitude independently, 
and does not depend on the relative numbers as a func¬ 
tion of magnitude. 


3.0.1. Causes of the Discrepancy 

A variety of tests were performed to identify the source 
of the discrepancy between the data and the available at¬ 
mosphere models at the faint end of the UV white dwarf 
cooling sequence. These included re-drizzling smaller 
sections of the image mosaic that excluded the dense 
core of 47 Tuc and performing PSF photometry on these 
regions, adjusting the sky subtraction algorithm to cor¬ 
rectly subtract the diffuse light between different frames, 
adjusting the annulus over which the sky is measured in 
the photometry to test for red leaks that were not be¬ 
ing accounted for, and exploring a range of apertures in 
a separate trial of aperture photometry analysis of the 
non crowded regions. We also performed a completely 
independent analysis of the data without the drizzling 
algorithm and with a different PSF-fitting photometric 
algorithm (J. Anderson, priv comm.). All of these tests 
resulted in colors for the faint white dwarfs consistent 
with our baseline reduction. 

Subsets of our core field mosaic overlap observations 
in other WFC3 UV filters, including F275W (GO-12311; 
PL G. Piotto) and F390W (GO-11664; PL T. Brown). 
We retrieved all of these data from MAST and subjected 
them to the same drizzling and PSF-fitting operations 
as our baseline data set, and matched the positions of 
sources across all four filters. While these data do not 
have the same depth and spatial coverage of our new ob¬ 
servations, they hint that the color offset with respect to 
the models may also appear in other color combinations 
that span a wavelength ~ 300nm (e.g., a possible offset 
is also seen in a F275W vs F3901U comparison). 

The source of this mismatch remains a mystery, and 
could be the consequence of one or more independent 
reasons. For example, the images in the redder filters 
like F3361U and F3901U have a very high level of crowd¬ 
ing as we approach the cluster center, and therefore the 
completeness is <75% even at 2 core radii. So, if there is 
a systematic issue with the photometry that is not being 
accounted for by the artificial star tests, then incomplete¬ 
ness could be the cause. Alternatively, it is possible that 
the observed colours are correct and rather the current 
generation of white dwarf atmosphere models lack an im¬ 
portant source of opacity at UV wavelengths. 

3.1. The Bolometric Luminosity Function 

The data presented here represent a new opportunity 
to construct a luminosity function for hot white dwarfs. 
This is particularly difficult in the field because the short 
cooling times imply that hot white dwarfs are relatively 
rare, requiring large scale surveys to characterise them. 
Observations in a globular cluster can isolate a well de¬ 
fined sample and better illuminate features of the lu¬ 
minosity function, at the price of foregoing some of the 
more detailed atmospheric characterisation possible from 
spectra of the closer field stars. In § [S] we will provide 
a detailed model fit to the data, but for the purposes 
of comparison to prior studies, it is also worthwhile to 
present the data in the same form as in other studies. 

The construction of an absolute bolometric magnitude 
requires the adoption of a cluster distance and redden¬ 
ing. Fortunately, the sparseness of the data in wide field 
surveys means we can compare results in relatively large 
(1 magnitude) bins, which renders the uncertainties in 
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the discussion above (of the order of 0.1 magnitudes) ir¬ 
relevant. We will adopt /tq = 13.3 and E{B — V)= 0.04, 
consistent with our determination above and our prior 
analyses of this cluster. In later sections, we will allow 
these quantities to vary in the fitting procedure, but we 
keep them fixed in this case because the uncertainties are 
small relative to our bin size. 

Our atmosphere model colours are based on pure Hy¬ 
drogen atmosphere models from Tremblay, Bergeron & 
Gianninas (2011), as a function of Te// and surface grav¬ 
ity. In the case of the WFC3 field, the F336W data 
are the deepest (photometric errors < 0.1 magnitude at 
a fainter position along the cooling sequence) and, in 
the case of the ACS fields, the F555W data are deepest. 
Thus, we convert the F336W and F555W magnitudes to 
bolometric magnitudes using the above models. The as¬ 
sumption of a pure Hydrogen atmosphere model for the 
hottest white dwarfs is questionable because it takes a 
finite time for the heavy elements to sediment out and 
for a star to take on the traditional onion-skin struc¬ 
ture assumed for white dwarfs. Fortunately, the spec¬ 
troscopic classifications of the SDSS hot white dwarfs 
from Krzesinski et al. (2009) clearly identify this tran¬ 
sition, showing that DA white dwarfs appear for bolo¬ 
metric magnitudes Mboi > 1.5. Thus, for magnitudes 
brighter than this, our bolometric magnitudes may be 
inaccurate, but the bulk of our sample is fainter than 
this. 

To construct the bolometric luminosity function in the 
outer fields, we must account for the effects of photomet¬ 
ric incompleteness. We first count and bin observed stars 
as a function of F555W to obtain a raw number count. 
In order to understand the incompleteness, we perform 
a monte carlo simulation in which we draw white dwarfs 
of the appropriate magnitude and then quantify their 
expected detectability using the results of the artificial 
star tests. This is sampled until the number of detec¬ 
tions equals the observed raw number, and the number 
of drawn input stars required to meet this target yields an 
estimate of the true underlying number of white dwarfs 
of this magnitude in our field. 

In the case of the luminosity function in the core, we 
must also account for the fact that the density of sources 
varies across the field, and this leads to a covariance be¬ 
tween incompleteness and location in the field. Our arti¬ 
ficial star tests determine incompleteness as a function of 
radius, and so we need to account for the radial profile of 
the underlying sources as well. Fortunately, 47 Tuc has 
been a subject of close study for many years, and the 
radial profile of stars is known to be well described by 
a King model with Wq = 8.6 (McLaughlin et al. 2006). 
Our core field extends out to 2.7', and so we are able 
to measure a core radius for MS stars directly from our 
data, yielding rg = 29.7" ± 0.4". This is within 2tT of 
that measured by Goldsbury, Heyl & Richer (2014) from 
surface density profiles. 

The core radii of white dwarfs and main sequence stars 
may not be the same, however, because of the mass loss 
associated with the latter stages of stellar evolution. The 
energy equipartition in such a dynamically relaxed sys¬ 
tem implies that lower mass stars such as white dwarfs 
(~ 0.5Mq) should have a larger velocity dispersion and 
therefore larger radius than more massive main sequence 
turnoff stars (~ 0.8Mq). Thus, our modelling needs to 


consider a range of possible core radii. Furthermore, 
there is a distinct possibility that the core radius may 
be a function of white dwarf luminosity. If the core re¬ 
laxation time is longer than the timescale over which a 
star loses mass, then the hottest white dwarfs may still 
exhibit the core radius of their more massive progeni¬ 
tor population. This results in an intimate connection 
between the details of stellar evolution and the number 
counts of hot white dwarfs. If much of the mass is lost 
on the Asymptotic Giant Branch, then the hottest white 
dwarfs may be over-represented in the core relative to 
their cooler brethren, as they will still exhibit a smaller 
velocity dispersion. Alternatively, if much of the mass is 
lost on the red giant branch, then we expect a smaller 
effect of this sort, because the stars will have spent ^ 10 ® 
years on the Horizontal branch with a mass close to that 
of a white dwarf, and so will have adjusted their ve¬ 
locity dispersion appropriately before ever becoming a 
white dwarf. Thus, we wish to examine a range of core 
radii ranging from one similar to the MS stars, up to 
^ 1.27 X 29.7" ^ 38", based on a white dwarf mass of 
O. 53 M 0 and an MSTO mass ^ O. 86 M 0 (e.g. Thompson 
et al. 2010 ). 

The construction of our resulting empirical luminos¬ 
ity function proceeds as follows when using the F336W 
data, mapped onto bolometric magnitudes using the at¬ 
mospheric models. In order to keep our results less sen¬ 
sitive to the modelling, we restrict our attention to radii 

> 30" from the core, in order to avoid uncertainties as¬ 
sociated with the aforementioned dynamical adjustment. 
On these scales, the two-body relaxation time locally is 

> 10 ® years, and therefore should not affect the lumi¬ 
nosity function at magnitudes appropriate to neutrino 
cooling. We count and bin observed stars as a func¬ 
tion of F336W to obtain a raw number count. In order 
to understand the incompleteness, we perform a monte 
carlo simulation in which we draw white dwarfs of the 
appropriate magnitude from the underlying King model, 
to form a radial profile and then quantify their expected 
detectability using the results of the artificial star tests. 
This is sampled until the number of detections equals 
the observed raw number, and the number of drawn in¬ 
put stars required to meet this target yields an estimate 
of the true underlying number of white dwarfs of this 
magnitude in our field. The error bars are calculated by 
repeating the procedure for a range of potential core radii 
from 20”-40” (to account for uncertainties in the cluster 
profile) and also by accounting for the Poisson fluctua¬ 
tions in the raw counts. In this manner, we define an 
empirical bolometric white dwarf luminosity function in 
the core of 47 Tucanae while accounting for photometric 
scatter and uncertainties in the mass profile. 

The resulting luminosity functions from the two data 
sets are given in Table [T] Figure [5] also compares these 
two independent determinations to the luminosity func¬ 
tions derived for the solar neighbourhood from the SDSS 
by De Gennaro et al. (2008) and Krzesinski et al. (2009). 
In those cases, a correction for the number density of 
stars has to be made, so we have allowed for an over¬ 
all normalisation shift and have set the normalisation so 
that the disk and cluster LF do not overlap exactly, to 

® The potential effects of dynamical relaxation are dealt with in 
more detail in Heyl et al 2015. 
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Fig. 5. — The filled circles indicate the bolometric function mea¬ 
sured for the inner field using the WFC3 data, and the open circles 
indicate the equivalent measured using the ACS data in the outer 
fields. The similarity in the number counts is somewhat coinci¬ 
dental, resulting from the competing effects of lower stellar density 
but greater field coverage in the ACS data. Also shown, as trian¬ 
gles, are the equivalent luminosity functions measured for the disk 
by SDSS survey, as presented by De Genarro et al. (2008) - open 
symbols - and Krzesinski et al. (2009) - filled symbols. These den¬ 
sities have each been multiplied by an overall normalisation factor 
to convert them to a number on this plot so that one can visually 
compare the resulting slopes. 

enable ease of visual comparison. The Krzesinki et al. 
(2009) luminosity function is the only one that goes to 
bright enough magnitudes to provide a meaningful com¬ 
parison with our data set. The precise number counts at 
the bright end of that data are somewhat uncertain be¬ 
cause of several completeness effects related to the spec¬ 
troscopic selection, as described in Krzesinski et al. and 
in Eisenstein et al. (2006). However, the ability to get 
spectroscopy of their stars allows them to classify stars 
according to spectral type, for which our stars are too 
faint. As noted above, genuine DA white dwarfs only 
appear at M^oi > 1.5 in the SDSS sample, as hotter 
white dwarfs have not yet had sufficient time for grav¬ 
itational settling to occur, and so our bolometric cor¬ 
rections are likely to be inaccurate for M^oi < 1. The 
slopes of the disk and cluster luminosity functions com¬ 
pare well for Mboi > 5 though, with the exception of a 
possibly larger flatening observed in the ACS data for 
Afboi > 8. At these magnitudes, the scatter in the white 
dwarf sequence and the SMC main sequence begins to 
overlap, and the lower counts are possibly the result of 
white dwarfs being lost due to confusion with SMC main 
sequence stars. 

At the bright end, our luminosity functions are consis¬ 
tent with a relatively smooth extrapolation of the slope 
at the fainter end towards higher luminosities. This con¬ 
tradicts the expectations based on the SDSS luminosity 
function. If we normalise the Krzesinski et al. LF to ours 
at Mboi = 4, the SDSS data would predict 10-19 stars 
(depending on which SDSS LF binning is used) in the 
magnitude bin centered on M^oi = 2, whereas we count 
only two stars in this bin, a discrepancy of approximately 
3 a. Torres et al. (2014) have commented on the SDSS 


TABLE 1 

Bolometric Luminosity Functions for 47 Tucanae. 




F336W 


F555W 

Mboi 

N 

N’ 

N 

N’ 

0 

1 

1±1 

0 


1 

1 

1±1 

0 


2 

2 

9 + 2 
^-2 

0 


3 

7 


0 


4 

8 

loli 

3 

0 + 2 
'J_2 

5 

20 

27+1 

12 

1311 

6 

41 

5S+11 

21 

23l® 

7 

89 


72 
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269 

417t|^ 

214 

25111° 

9 

513 

845t^^ 

279 

350111 

10 

809 

1436):^(® 

405 

551111 


Note. — Under each bandpass we list the raw counts (N) in 
each bolometric magnitude bin and then the counts corrected for 
selection biases {N') as described in the text. The error bars in¬ 
corporate the effects of photometric scatter, incompleteness and 
the radial density profile, and are thus moderately larger than the 
nominal Poisson error. Note that the bolometric corrections for 
the brightest two bins (Af(,oi < 1.5) may be uncertain as these may 
be too young to be pure H atmospheres. 

‘plateau’ and found that it could not be reproduced with 
standard models. Torres et al examined a variety of pos¬ 
sible explanations and concluded that the upturn cannot 
be the result of a feature in the Galactic inputs such as 
star formation rate variations. They further concluded 
that errors in the model colours, uncertain white dwarf 
cooling tracks or statistical fluctuations were unlikely to 
be the cause. Their preferred solution was that the ob¬ 
served flatening was the result of erroneous mass deter¬ 
minations that led to some white dwarfs being assigned 
unrealistically low masses and consequently larger dis¬ 
tances and larger inferred luminosities. All of our white 
dwarfs are at the same distance and have the same mass, 
since they follow a single cooling track, so that this selec¬ 
tion effect should not influence this luminosity function. 
The uncertainty in the bolometric corrections for the core 
data should also not affect this conclusion, because the 
bins are quite large, and any white dwarfs observed would 
be counted. In summary, the lack of any plateau in our 
data confirms the conclusion of Torres et al. that this is 
likely the result of systematic error rather than a phys¬ 
ical effect. This is consistent with the behaviour of the 
proper motion selected LF of Rowell & Hambly (2011) 
at the bright end (shown in Figure 3 of Rowell 2013), 
which also drops monotonically except for a spike at the 
brightest magnitudes which may represent the presence 
of non-DA white dwarfs. 

4. MODELS 

In order to evaluate the role of various physical inputs 
into the cooling of white dwarfs, we need to match our 
observed luminosity function to models for the rate of 
cooling. We have established an extensive framework for 
such modelling in previous papers on the cooling func¬ 
tions of cluster white dwarfs, most recently Hansen et 
al. (2007) and Hansen et al. (2013). However, our previ¬ 
ous work has been devoted primarily to the coolest white 
dwarfs and the focus on hot white dwarfs in this project 
necessitates the investigation of several relevant factors 
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(g.cm 3) 

Fig. 6.— The solid curve indicates the evolution of the cen¬ 
tral density and temperature of a 0.9Mq star with standard input 
physics, while the dashed line indicates the equivalent but for a 
star without neutrino cooling. The two large solid points indicate 
the location at which each star passes through an effective temper¬ 
ature of lO^K - an approximate start of the white dwarf cooling 
sequence. The dotted contours indicate lines of constant neutrino 
emissivity, increasing by two orders of magnitude with each contour 
as one moves from the lower left to the upper right. 

with greater care and detail than in our previous work. 

Our default white dwarf cooling models are based on 
the code of Hansen (1999). For the purposes of this work, 
the most important part is the implementation of neu¬ 
trino cooling, which is taken from the work of Itoh et 
al. (1996). The focus on young, hot, white dwarfs also 
requires greater attention to the starting conditions of 
stars that begin to descend the white dwarf cooling se¬ 
quence. In order to establish a baseline starting model, 
we have used the MESA models (Paxton et al. 2011) 
to construct full evolutionary models of stars with initial 
masses ^ 0.9Mq and metallicity 0.004. Figure [S] shows 
the evolution of the central density and temperature for 
two stellar models of initial mass 0.9Mq and metallic¬ 
ity 0.004, both with (solid curve) and without (dashed 
curve) neutrino emission. The dotted contours indicate 
loci of constant neutrino emissivity. The two solid points 
indicate the location at which the star reaches an effec¬ 
tive temperature of 10® K, an approximate location for 
the start of the white dwarf cooling sequence. These 
curves indicate that a self-consistent treatment of neu¬ 
trino cooling in the stellar evolution implies differences 
of a factor of two in the initial white dwarf central tem¬ 
perature. They also indicate that the dominant neutrino 
emission process of relevance here is plasmon emission. 
The change in slope of the dotted contours corresponds to 
the transition from emission due to photoneutrino pro¬ 
duction (left of the diagram) to plasmon decay (lower 
right). The white dwarfs are well into this latter regime. 

We use these models as the reference model and con¬ 
sider variations in the input parameters in our cooling 
models to investigate the sensitivity to uncertain aspects 
of the physics. 

The most obvious parameter of relevance is the total 
white dwarf mass. We consider a range of masses ^0.5- 


0.55Mq. This is the mass range expected on theoretical 
grounds (e.g. Fusi-Pecci & Renzini 1976) and observed 
directly from the pressure broadening of spectra of hot 
cluster white dwarfs (Moehler et al. 2004; Kalirai et al. 
2009). Modelling of the full white dwarf cooling sequence 
in this cluster and others (Hansen et al. 2007; 2013) also 
agrees with this estimate. In that case, the mass can 
be constrained directly because the observations contain 
mass information due to the fact that the onset of crys¬ 
tallisation is mass-dependant and so the cooling curves 
for different masses are not scaled versions of one an¬ 
other. The sample of white dwarfs observed here is too 
hot to have undergone crystallisation and so we restrict 
our analysis to the plausible mass range above identified 
by other means. 

In addition, although we start with a reference model, 
there is also some uncertainty in the establishment of 
this baseline. The chemical composition of the white 
dwarf core is affected by prior evolution as the temper¬ 
ature history of the core, as well as uncertainties in the 
cross-section of some nuclear reactions, will affect the rel¬ 
ative abundances of Carbon and Oxygen and their profile 
through the core. We have studied this in detail using 
a range of stellar models, as detailed in the supplemen¬ 
tary information of Hansen et al. (2013). These models 
also account for possible variations in the initial Helium 
abundance of the star and uncertainties related to possi¬ 
ble mixing of compositional gradients by Ledoux convec¬ 
tion. Incorporating these uncertainties produces a range 
in central Oxygen mass fraction of 0.64-0.86. The core 
composition proves to be much less important for the 
hot white dwarfs discussed here than for the cooler white 
dwarfs undergoing crystallisation, because the neutrino 
emission is quite insensitive to the relative proportions of 
Carbon and Oxygen in the core. Furthermore, although 
the MESA models provide an initial temperature profile, 
with a central temperature ~ 10®K with neutrino cooling 
and ~ 2 X 10®K without, we also allow for a variation of 
10% in starting temperature to account for uncertainties 
in prior evolutionary stages. 

We have also examined a range of Helium and Hydro¬ 
gen layer masses. These lie at sufficiently low densities 
that they will not have an effect on the neutrino emis¬ 
sion, but can contribute to the regulation of thermal dif¬ 
fusion of energy through the white dwarf and can there¬ 
fore potentially affect the location of the transition from 
neutrino-dominated to photon-dominated cooling. Hy¬ 
drogen layers can have two additional effects. The lower 
mean molecular weight of Hydrogen means that the ra¬ 
dius of the white dwarf can be non-negligibly affected by 
the thick hydrogen layer, at least for hot white dwarfs. 
This can change the location of the cooling sequence in 
the colour-magnitude diagram and therefore affect the 
fit of the data. Furthermore, sufficiently thick Hydrogen 
layers can result in a contribution to the luminosity from 
p-p nuclear burning, which can also change the cooling 
curve. 

For young white dwarfs, it is also possible that the H 
and He layers have not completely gravitationally sepa¬ 
rated, which will change the opacity of the atmosphere 
and affect the cooling. This occurs near the surface of the 
white dwarf and will therefore have a negligible effect on 
the neutrino cooling, but it may affect the location where 
photon cooling dominates over neutrinos. We note, how- 
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Fig. 7.— The solid curve shows the photon luminosity of a 
O.53M0 white dwarf with standard neutrino cooling. The corre¬ 
sponding model with no neutrino luminosity is given by the dotted 
curve. The dashed curve represents the same model except that 
we have increased the initial H layer mass. 

ever, that we examine a large range of Hydrogen and 
Helium layer masses, and that the Helium layer masses 
are an order of magnitude or more larger than the Hy¬ 
drogen layer masses. Thus, a surface layer of Hydrogen 
diluted with Helium can be considered the limiting case 
of a thin Hydrogen atmosphere. We have considered H 
layer mass fractions as small as qn = 10“® and find that 
our results are very insensitive to values of qh < lO”"*. 
Thus, we believe that uncertainties in the separation of 
H and He do not affect our results. 

4.1. Binary White Dwarfs 

The WFC3 data are taken in the cluster core and so 
any effects of binarity are likely to be increased due to 
mass segregation. However, the hot white dwarf sequence 
observed here covers a very short-lived phase of stellar 
evolution and so the observation of two hot white dwarfs 
in the same binary would require that the original bi¬ 
nary pair had main sequence masses within 1% of each 
other. If the mass ratio is determined by randomly sam¬ 
pling the Salpeter mass function, this occurs in ~ 10“^ of 
pairs containing one main sequence turnoff star. Thus, 
even for relatively high global binary fractions, we do 
not expect a substantial contribution of binaries to this 
luminosity function. 

4.2. Non-standard neutrino physics 

There are many proposals for how modifications to the 
properties of the neutrino could point the way to physics 
beyond that contained in the standard model. One of 
the most prominent is the proposal that the neutrino 
may possess a non-zero magnetic dipole moment. Such 
a modification also has an effect on the plasmon neu¬ 
trino cooling rate and is thus amenable to constraint in 
hot white dwarfs (e.g. Wang 1992, Blinnikov & Dunina- 
Barkovskaya 1994, Haft, Raffelt & Weiss 1994). 

We investigate the resulting extra cooling by consider¬ 
ing the inclusion of an amplification factor for the plas¬ 


mon neutrino emission from Haft et al. (1994) 


^total 

^standard 


l + 0.318g 


lO-iVs 


/ hujp \ 

\10keV J 


(1) 


where jii, is expressed in terms of the Bohr magneton /i b 
and LOp is the plasma frequency, given by 


hujp = 0.021keVp^^^ 1 + (p/1.96 x lO^g/cm 

( 2 ) 

Q is a function that depends on the underlying plasma 
properties but is ^ 1 to within 10% over the density and 
temperature range appropriate to the centers of white 
dwarfs. 


4.3. Axion Cooling 

Hot white dwarfs are a favourite astrophysical con¬ 
straint on hypotheses of exotic particles, because such 
particles can cause extra cooling if they couple to nu¬ 
clei, electrons or photons. One commonly proposed extra 
cooling mechanism is emission of axions. The axion is a 
hypothetical particle (Weinberg 1978; Wilczek 1978) as¬ 
sociated with the U(l) symmetry proposed by Peccei & 
Quinn (1977) as a solution to the strong CP problem of 
particle physics. If such axions couple to electrons they 
can result in an increase in the luminosity of hot white 
dwarfs. 

We can treat this by including an extra cooling contri¬ 
bution in the white dwarf models, of the form 

ea. = 1.08 X 1023ergsg-is-iF(p,T)^^r/, (3) 

where Ty is the temperature measured in units of lO^K, 
and F is the correction for ionic correlation effects from 
Nakagawa, Kohyama & Itoh (1987). The factor Z'^/A 
ranges from 3.6 to 3.8, depending on the core composi¬ 
tion. The axion-electron coupling gae is related to the 
mass of the axion max by gae = 2.8 x \Q~^‘^max/^'meV 
(up to an unknown angle (3). 

4.4. Radiation into Extra Dimensions 

Some models of gravitation invoke a version of string 
theory that includes additional dimensions that are com- 
pactified on dimensions large enough to have experimen¬ 
tal consequences (Arkani-Hamed, Dimopolos & Dvali 
1998). In such models, the coupling between photons, 
electrons or nucleons to gravity can result in the emis¬ 
sion of Kaluza-Klein gravitons into these extra dimen¬ 
sions, and will manifest themselves as an extra source 
of cooling for hot stars. Barger et al (1999) present ex¬ 
pressions for the emissivity for a variety of processes and 
place constraints on the mass scale Mg corresponding 
to extra dimensions by using astrophysical constraints 
from the Sun, red giants, and supernova 1987A. Biesi- 
ada & Malec (2002) place additional constraints on this 
mass scale using the pulsational stability of the pulsat¬ 
ing white G117-B15A, and the fact that electron-graviton 
bremstrahlung would accelerate the cooling. 

We incorporate the emissivities of Barger et al for three 
processes, namely photon-photon annihilation, Gravi- 
Compton-Primakoff scattering, and gravibremstrahlung. 
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These are 

= 5.1 X 10-®ergs p-^ 

eGCP = 4.5 X 10"®ergsg"^s“^T/ (5) 

eGB = 5.8 X lO-Vgsg-is-'Z^T^ , (6) 

where pg is the density in units of 10 ®g.cm“^, Ty is the 
temperature in units of 10^ K and Zj is the mean ion 
charge relative to Nitrogen (the value for our model mix¬ 
tures of Carbon and Oxygen range from 7.2 to 7.8). 
We have assumed two macroscopic extra dimension^. 
These rates are calculated for non-degenerate material, 
although Barger et al. claim that corrections for degener¬ 
acy should be of order unity. In the case of the red giant 
cores considered in that paper, such an approximation 
may be acceptable, but it is likely to be an overestimate 
under the much more degenerate conditions appropriate 
to a white dwarf (a fact that does not appear to have 
been acknowledged by Biesiada & Malec). The latter 
two processes above depend on electron scattering and 
therefore the rates should be considerably depressed in a 
white dwarf. Although a full rederivation of the rates is 
beyond this paper, we also test modified expressions for 
CGCP and eoB in which the rates are scaled by a factor 
kT/Ep, where Ep is the Fermi energy of the electrons. 
This limits the fraction of the electrons that participate 
to only those within kT of the Fermi surface, and which 
therefore are energetically able to scatter to a state that 
is not excluded by the Pauli principle. The resulting ex¬ 
pressions are 

CGCP = 2.3 X lO-^ergs g-^s-^T^pf'^ (l^) 

CGB = 3 X 10“®ergs g-^-^Z'^T^ (l^^) 

4.5. Monte Carlo Model 

Using the cooling curves constructed in this fashion, 
we simulate artificial white dwarf populations assuming 
that the birth rate of white dwarfs has been constant 
over the last 10 ® years (since our comparison data set 
comprises ages ^ 10® years at the faint end). For each 
white dwarf, the model provides a unique luminosity and 
effective temperature, which is converted into a magni¬ 
tude and colour based on a chosen value of distance and 
extinction. We then use the results of our radial density 
profiles and photometric scattering matrix to probabilis¬ 
tically assign output magnitudes and colours given the 
original input value, including the nondetection of stars. 

The resulting model populations are then binned in 
colour and magnitude on the grid shown in FigurelU and 
compared to the data binned in the same manner. Given 
the problematic colours for the core data, we restrict this 
analysis to the swathe of outer cluster fields observed 

^ One extra dimension is excluded experimentally by confirma¬ 
tion of Newtons gravitational force law and three dimensions pro¬ 
duce rates too low to be tested here. 
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Fig. 8. — The points show the observed white dwarf cooling se¬ 
quence in the ACS bandpass, while the red boxes indicate the grid 
used to bin the data and models for the purposes of a quantitative 
comparison. 

with the ACS camera. The grid is chosen to provide 
sufficient colour resolution to constrain the distance and 
extinction, while still maintaining enough points to en¬ 
sure stability of the y® fitting procedure. The brightest 
bin covers 1 magnitude in range because the rapid cool¬ 
ing at these temperatures makes the counts insensitive 
to the cooling and more sensitive to initial conditions. 

We restrict our attention to the stars with F555W< 
25.15. Our limit is set by the increasing risk of con¬ 
fusion due to photometric scatter from the SMC main 
sequence in the background. We also account for this in 
our reddest bins by modelling the colour spread of the 
SMC population as a function of magnitude. We extrap¬ 
olate this model bluewards with a maxwellian function, 
using the fitted colour dispersion, in order to anticipate 
the pollution of the white dwarf sequence with scattered 
main sequence stars, and include these in the model. 

5. RESULTS 

We compare our data with models calculated as de¬ 
scribed in § m using the reference model from MESA 
and considering variations in initial temperature and core 
chemical composition. Eor each core model, the com¬ 
parison with the data also requires a choice of distance 
modulus and extinction. These, combined with the pa¬ 
rameters of white dwarf mass, and Helium and Hydrogen 
surface mass fractions, as well as an overall global nor¬ 
malisation, imply a total of six parameters to be varied 
in the comparison of a particular model with the data. 
In all cases, we allow the distance modulus to vary from 
13.0 to 13.5 and the reddening to vary from E(B-V)=0.0 
to 0.08. White dwarf masses are varied from 0.5Mq to 
O. 6 A/ 0 , while Helium and Hydrogen mass fractions are 
varied from qpe = 0.01-0.04 and = 1 x 10“®-6 x 10“^. 

The comparison between models and data is performed 
by binning on the grid shown in Figure [5] The brightest 
stars are assigned to a single large bin in order to lessen 
the sensitivity of our results to uncertainties in initial 
parameters. The principal goal of this comparison is to 
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Fig. 9.— The points represent the observed counts in each bin 
of the grid shown in Figure [8] The points are ordered by colour 
from blue to red, within a given magnitude bin, and the magnitude 
increases from left to right. The histogram (shown in red) is the 
best fit model described in the text. 

locate the magnitude at which the cooling rate changes 
due to the transition to photon-dominated cooling from 
other processes, be it via neutrinos or more exotic pro¬ 
cesses. The colour binning is chosen to constrain the 
location of the cooling sequence in colour and thereby 
measure the extinction simultaneously with the cooling. 
The final grid thus contains 47 bins, so that comparison 
with a particular core model has 40 degrees of freedom. 
The best fit for a model with our standard neutrino cool¬ 
ing yields = 44.8, for O.53M0, a Helium layer mass 
fraction qne = 0.037 and ( 7 _y = 4 x 10“^. This Hydro¬ 
gen layer is slightly larger than usually assumed, but not 
so large that residual Hydrogen burning is a major con¬ 
tributor. Larger masses {qu > 6 x 10“^) yield poor fits 
(x^ > 100) because the nuclear contribution dramati¬ 
cally changes the shape of the cooling curve, indicating 
that the cooling of the white dwarfs in 47 Tuc are not 
being substantially delayed by residual H burning, as has 
been suggested for lower metallicity progenitors (Miller 
Bertolami, Althaus & Garcia-Berro 2013). The best fit 
distance modulus is no = 13.27, which compares well 
with other distance estimates, although the best fit red¬ 
dening, E(B-V)=0.08, is at the red end of the allowed 
range. Figure IHl shows the comparison of the model with 
the data. In this plot, we show the counts as a function of 
colour in each luminosity bin (separated by vertical dot¬ 
ted lines). For instance, the five bins between the dotted 
lines indicating F555W=24.1 and 24.25 show the counts 
in five colour bins spanning F4351F — F5551F = —0.1 to 
0.15. This projection of the two-dimensional data into a 
pseudo-luminosity function allows us to assess the quality 
of the model fit to the data. 

We test the sensitivity of our results to the neutrino 
emission by applying a uniform scaling factor fs to 
the neutrino emission and redoing the fit, marginalis¬ 
ing over the white dwarf mass, internal composition, and 
layer masses, as well as distance modulus and extinction 
ranges described above. The resulting curve of x^ as a 
function of fs is shown in Figure [TOl At 95% confidence. 
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Fig. 10.— The solid points show the minimum each value 

of /s, after marginalising over white dwarf masses and values of 
q^e QHi well as distance and extinction. The dotted line 
shows the value of Ax^ = 4 relative to the minimum, so that the 
dashed vertical lines indicate the region of 95% confidence. The 
slope is steeper for fs < 1 than for /s > 1, suggesting that the 
constraints on a modicum of extra cooling are weaker than models 
with less cooling. 

0.6 < fs < 1.7, with the best fit x^ achieved for fs ^ 0.9. 
The value at /^ = 1 is only Ax^ = 0.5 larger than the 
minimum, so that the difference is not statistically sig¬ 
nificant. Thus, our white dwarf cooling models favour 
the traditional implementation of neutrino emission and 
constrain variations to be approximately within a factor 
of two. There is also no indication that increased neu¬ 
trino cooling can be offset by nuclear burning in a thick 
H envelope. 

The fs scaling models all start with the initial condi¬ 
tions determined from the MFSA models, but we noted 
in § m that a self-consistent evolutionary model with no 
neutrinos would begin with a higher central temperature, 
so we have also tested a model with no neutrino cooling 
starting with the appropriate initial conditions. The best 
fit is achieved for a O.55M0 model with qne = 0.025 and 
qn = 2 X 10 “"^. The best fit distance and extinction is 
fio = 13.16 and F(B-V)=0.08. Nevertheless, the mini¬ 
mum x^ = 184, which is very poor compared to our best 
fit model. Figure [Tl] shows this model and demonstrates, 
when compared with Figure [51 how this data discrimi¬ 
nates between models. In this case, the lack of neutrino 
cooling results in slower early cooling and an excess of 
counts at bright magnitudes relative to faint magnitudes. 

These results also allow us to anticipate a comparison 
with constraints on neutrino emission from the pulsa- 
tional stability of hot DBV white dwarfs (Winget et al. 
2004). Bischoff-Kim (2008) provides a model for the pul- 
sational properties of the hot DBV star EC20058-5234, 
and predicts the rate of change of two pulsational periods 
as a function of A (scaled explicitly relative to the plas- 
mon neutrino rate in that case, but effectively the same 
as our overall scaling parameter fs in this regime). Based 
on our constraints above and the figures in Bischoff- 
Kim (2008), we anticipate that the period derivatives 
of the 257s and 281s modes in that star should be within 
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Fig. 11. — The points represent the same observed counts as in 
Figure 0 but the red histogram is now the best fit model for the 
case of no neutrino cooling. We see that the requirement that the 
model fit the data at the faint end results in a large overestimate 
of the number of bright white dwarfs. 

1.0±0.4 X 10“^^s/s and 1.4±0.5 x 10“^^s/s respectively. 

5.1. The Neutrino Magnetic Moment 

The restriction on scaled neutrino emissivities should 
then also translate into constraints on exotic radiation 
mechanisms. The constraint fg < 1.7 can be used 
in equation 0 to obtain an approximate limit < 

5 X 10~^^fiB- A more accurate limit is obtained by re¬ 
peating the cooling calculations using the correction to 
the emissivity given by equation o, including the den¬ 
sity dependance. If we repeat the fitting procedure with 
these models, marginalising again over distance, extinc¬ 
tion and white dwarf mass parameters, we find that the 
95% confidence limit on the neutrino magnetic moment 
is /Ti/ < 3.4 X 

This limit is an order of magnitude better than the cur¬ 
rent laboratory limits ~ 3 x 10~^^/j,b (Daraktcheva et al. 
2005; Arpesella et al. 2008; Beda et al. 2010) and also 
better than astrophysical limits based on the lifetime of 
the Sun (Bernstein et al. 1963; Raffelt 1999) and Hori¬ 
zontal Branch Stars (Raffelt et al. 1989). Most directly, 
it represents a factor of several improvement of earlier 
determinations /ii 2 < 10 based on the hot white dwarf 
luminosity function (Blinnikov & Dunina-Barkovskaya 
1994). Although measurements of the hot white dwarf 
luminosity function have improved since then, the excess 
counts at the bright end (see discussion in S IXTI) have 
hampered the statistical significance of such constraints 
using the more up-to-date data (Miller Bertolami 2014). 
Our measurement is also comparable to the limits de¬ 
rived from the core mass of red giant stars (Castellani 

6 degl’Innocenti 1993; Catelan et al. 1996; Haft et al. 
1994; Raffelt 1990; Raffelt & Weiss 1992; Viaux et al. 
2013) but has the advantage of being measured in stars 
whose properties are dominated by the neutrino emission 
and in a manner in which the distance (the largest source 
of error in the red giant measurements) is constrained si¬ 
multaneously. 

In addition to the consequences for particle physics, 



Fig. 12. — The points represent the same observed counts as in 
Figure |9] but the red histogram is now the best fit model for the 
case of axion cooling with an axion mass of 17 meV. Fitting the 
model at the faint end results in an underproduction of stars at 
brighter magnitudes. 

a neutrino magnetic moment in the range /ii 2 20-50 
could lead to qualitatively different late-time evolution 
for stars in the initial mass range ~9-17Mq (Heger et 
al. 2009). However, our limit here appears to preclude 
this eventuality. 

5.2. The Axion Mass 

In a similar fashion, models with additional cooling 
due to axion emission, given by equation ©, can be 
compared to the data to obtain a 95% conhdence limit 
of gae < 8.4 x 10“^'*, which implies an axion mass 
< 4.1mey. 

This constraint contradicts the claim of Corsico et al. 
(2012), whose modelling of the pulsating white dwarf 
G117-B15.A implied that the default rate of cooling 
was too slow to match the observed rate of change of 
the pulsation period. They invoked an extra cooling 
due to axion emission that implied an axion mass of 
max = 17.4:'^2 ymeV, a 5.8 cr discrepancy with respect 
to our limit. However, the Corsico result was based 
on the rate of change of a particular oscillation mode 
which they claimed was trapped in the Hydrogen enve¬ 
lope. The trapping strongly slows the expected rate of 
period change and suggests that future improvements in 
the asteroseismological modelling with untrapped modes 
may better match the observations. 

Figure [12] shows the best fit solution for the case of a 
17 meV neutrino (with a distance hq = 13.21 and E(B- 
V)=0.08). A comparison with Figure 0 demonstrates 
how too much cooling adversely affects the model ht - 
namely that the counts at brighter magnitudes are now 
too low with respect to fainter magnitudes (the opposite 
trend seen in Figure fTTl) . 

Our limits are not consistent at the 95% confidence 
level with the proposal that an axion mass ~ bmeV 
is motivated by the overall behaviour of the hot white 
dwarf luminosity function (Isern et al. 2008). Caution 
has been expressed about this result by Miller Bertolami 
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et al. (2014), who note that the strength of this conclu¬ 
sion varies depending on which observational luminosity 
function is used, and is further muddied by the aforemen¬ 
tioned uncertainties in the bright luminosity function. 
Our results do not require such an axion, as the inclu¬ 
sion of axion cooling only worsens the fit of the model, 
but are not strong enough to conclusively rule it out. 


5.3. Radiation into Extra Dimensions 

Testing models which include radiation into extra di¬ 
mensions place a 95% limit on the mass scale Mg > 
0.54TeH/c^ when we include our ad-hoc corrections for 
degeneracy (ignoring the degeneracy correction increases 
this to 0.9TeVjc^). The limits are different because the 
steep temperature dependance of means that this 
rate is irrelevant even for our hot white dwarfs, and only 
becomes important at energies corresponding to super¬ 
nova cores, so that the cooling of the white dwarfs is 
indeed regulated by the processes that depend on elec¬ 
tron scattering. These expressions are comparable to 
constraints based on the properties of the Sun but weaker 
than those based on red giants and supernovae (Barger 
et al. 1999, Cassisi et al. 2000). They are also weaker 
than the constraint derived by Biesiada & Malec (2002), 
although that estimate is likely overstated because the 
rates were not corrected for the effects of degeneracy. 


6. CONCLUSION 

Our results suggest that the cooling of hot white dwarfs 
is presently well understood. The bolometric luminosity 
function we infer from our hottest stars shows a mono¬ 
tonic decrease to the highest effective temperatures ob¬ 
served, and does not display the plateau observed in 
some prior studies. This supports recent claims that this 
feature was the result of systematic uncertainties rather 
than the consequence of unanticipated physical effects. 
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However, our observations in the core of the cluster 
result in a problematic match between the data and the 
models, in the sense that the model colours appear too 
blue relative to the data at fainter magnitudes. This sug¬ 
gests that either the photometric calibrations are system¬ 
atically biased by the crowding of the field or that the 
short wavelength opacities are not fully understood for 
white dwarfs of effective temperatures < 19,000 K. 

The atmosphere models do fit the data for the outer 
cluster fields, observed in longer wavelength bandpasses. 
We ht detailed cooling models to the number counts in 
these fields, which provide an excellent match to the ob¬ 
servations without requiring any substantial modifica¬ 
tions to the physics. Indeed, the comparison suggests 
that the neutrino emission models of Itoh et al. (1996) 
are accurate to within a factor of two (at 95% confi¬ 
dence). This match can also be used to constrain a 
variety of proposals for additional exotic cooling mech¬ 
anisms, and provides constraints on the existence of a 
neutrino magnetic moment, emission of axions, and pos¬ 
sible radiation of energy into extra dimensions. In each 
of these cases, the inclusion of extra cooling makes the 
Ht to the data worse, and provides little evidence for any 
exotic physics. 
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